Can the Jet steepen the hght curves of GRB afterglows? 
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ABSTRACT 

Beaming of relativistic ejecta in GRBs has been postulated by many authors 
in order to reduce the total GRB energy, thus it is very important to look for 
the observational evidence of beaming. Rhoads (1999) has pointed out that the 
dynamics of the blast wave, which is formed when the beamed ejecta sweeping 
the external medium, will be significantly modified by the sideways expansion 
due to the increased swept up matter. He claimed that shortly after the bulk 
Lorentz factor (F) of the blast wave drops below the inverse of the initial 
opening angle (^o) of the beamed ejecta, there will be a sharp break in the 
afterglow light curves. However, some other authors have performed numerical 
calculations and shown that the break of the light curve is weaker and much 
smoother than the one analytically predicted. In this paper we reanalyse the 
dynamical evolution of the jet blast wave, calculate the jet emission analytically, 
we find that the sharp break predicted by Rhoads will actually not exist, and 
for most cases the afterglow hght curve will almost not be affected by sideways 
expansion unless the beaming angle is extremely small. Wc demonstrate that 
only when 9q < 0.1, the afterglow light curves may be steepened by sideways 
expansion, and in fact there cannot be two breaks as claimed before. We have 
also constructed a simple numerical code to verify our conclusion. 

Subject headings: gamma-rays: bursts 
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1. Introduction 

The BeppoSAX results have revolutionized our understanding of GRBs by opening a 
window on X-ray, optical and radio afterglows (e.g. Costa et al. 1997; Piro et al. 1998; Van 
Paradijs et al. 1997). The decaying power- law long- wavelength afterglows are explained as 
the emission from a relativistic blast wave which decelerates when sweeping up interstellar 
medium. The dynamical evolution of GRB fireballs and the emission features have been 
studied by many authors (e.g. Sari 1997; Meszaros, Rees & Wijcrs 1998; Wei & Lu 1998a,b; 
Sari, Piran & Narayan 1998), most of them considered the fireball being isotropic. 

The discovery of GRB afterglow shows that GRBs are at cosmological distances. If 
so, the total energy for typical GRB event is about 10^^ ergs. This year, an extraordinary 
event GRB990123 was detected, which was the brightest burst ever detected by BeppoSAX 
satellite, and is in the top 0.3% of all bursts (Feroci et al. 1999). The detection of the 
redshift showed that the burst appears at ^ > 1.6, with its 7-ray fluence of ~ 5 x 10~^ 
ergs cm~^, the total energy of this source is > 1.6 x 10^*^ ergs if the emission is isotropic 
(Andersen et al. 1999; Kulkarni et al. 1999). This energy is so large that it gives a great 
challenge to the popular models. For models involving stellar mass central engines it is 
necessary to assume that the ejecta are beamed in order to explain such a huge energy. 

Now the main uncertainty of bursts' energy is whether the bursts radiate isotropically 
or are beamed into a small solid angle. As shown above, the extreme large energy favors 
the emission being beaming. Then there is one question: how can we identify the radiation 
between beaming or not? Rhoads (1997, 1999) have shown that the lateral expansion of 
the shocked, relativistic plasma causes that at some moment the surface of the blast wave 
starts to increase faster than due to the cone-outflow alone, then the blast wave begins to 
decelerate faster than without the sideways expansion since more interstellar medium has 
been swept up by blast wave. Rhoads claimed that this effect will produce a sharp break in 
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the GRB afterglow light curves. Such a break is claimed to be present in the light curves of 
GRB990123 and GRB990510 (Kulkarni et al. 1999; Harrison et al. 1999). Sari et al. (1999) 
speculate that afterglows with very steep light curves are highly beamed. 

However, some other authors (Panaitescu & Meszaros 1998; Moderski et al. 1999) have 
performed numerical calculations of dynamical evolution of blast wave, and shown that the 
break of the light curve is weaker and much smoother than the one analytically predicted. 
Thus there are two opposite conclusions about the jet emission, the analytical treatment 
predicts a sharp break, while the numerical calculation shows no such sharp break. 

In this paper we will first give an analytical treatment of the dynamical evolution of 
the jet blast wave and its emission features, and will demonstrate that the sharp break will 
not actually exist, we may observe the steepening of the light curve only when the jet angle 
is extremely small, i.e. < 0.1. We then perform a simple numerical calculation to confirm 
our results. In next section we discuss the dynamical evolution of blast wave, in section 3 we 
calculate the jet emission analytically, and finally we give some discussions and conclusions. 



2. Dynamical evolution of the jet 

Now we consider an adiabatic relativistic jet expanding in surrounding medium. For 
energy conservation, the evolution equation is 

rV = const (1) 

where F is the bulk Lorentz factor, and V is the jet volume, V — 27rr^(l — cos9j)/3 oc r^Oj 
for 9j <^ 1, and 9j — 9o + 9' — Oq + Cgtco/ct, where 9o is the initial jet opening half-angle, 
9' describes the lateral expansion, Cg is the expanding velocity of ejecta material in its 
comoving frame, and t (tco) is the time measured in the burster frame (comoving frame). For 
relativistic expanding material it is appropriate to take Cg to be the sound speed Cg — 
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(Rhoads 1997, 1999), and Rhoads (1999) has given tco/t = 2/bV. Since the jet expands 
relativistically, there is the relation T oc r/F^, where T is the time measured in the observer 
frame, r = ct is the radial coordinate in the burster frame. From above relations, we have 

r(i + Ej)V4 ^ T-'/s (2) 

where F^ = —Oq-^. In this paper all quantities with the subscript "b" denotes the point 
at which Cstco = r9o, which means that after that time the sideways expansion begins to 
dominate the radial divergence. Since T oc r/F^, so F;,/F = {^^V^^- For the case T ^ Tj,, 
it is well known that F oc T'^/^, then r oc F^T oc T^/^, and F^/F = {T/Tkf/^. For another 
case T 3> T^, the radial coordinate r is nearly a constant, r ~ r^, so F^/F ~ (T/T^)^/^. 
Therefore we have 

\ T-3/8[l + (^)V2]-i/4^ iiT>n 
It is obviously that F oc T-^/^ for T <C Tt, and F oc T'^/^ for T > Tf,. The rapid decrease 
with time of Lorentz factor F is due to the fact that larger amounts of surrounding matter 
has been swept up by ejecta (Rhoads 1997, 1999). 

In the following we calculate the value of Tf,. According to the fireball model, the 
decelerating radius of the ejecta is 

" ^Tr^gFgmmpC^^'^' 

where E is the burst energy, and Fq = E/Mqc^, Mq is the initial baryon mass. Rhoads 
(1999) has given ri, — [-^^^j^Y^^i^d- The relation between T and r is T = r/^F^c, where 
the numerical value of C lies between ~ 3 and ~ 7 depending on the details of the 
hydrodynamical evolution and the spectrum. Sari (1997) and Waxman (1997) have shown 
that the typical value of ( is about 4, then the time 
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We see that the break time Tj, is very large for typical parameters, which means that the 
transition from F oc y-^/s F oc T~^/^ is usually very slowly and smoothly 



Now we calculate the emission flux from the jet. Here we adopt the formulation and 
notations of Mao & YI (1994). In our model the ejecta is flowing outwards relativistically 
(with Lorentz factor F) in a cone with opening half angle 9j. For simplicity, we assume that 
the radiation is isotropic in the comoving frame of the ejecta and has no dependence on the 
angular positions within the cone. The radiation cone is uniquely defined by the angular 
spherical coordinates (^,0) of its symmetry axis, here 9 is the angle between the line of sight 
(along z-axis) and the symmetry axis, and is the azimuthal angle. Because of cylindrical 
symmetry, we can assume that the symmetry axis of the cone is in the y — z plane. In 
order to see more clearly, let us establish an auxiliary coordinate system {x', y', z') with the 
2;'-axis along the symmetry axis of the cone and the x' parallel the x-axis. Then the position 
within the cone is specified by its angular spherical coordinates 9' and (0 < ^' < 9j, 
< 0' < 2tt). It can be shown that the angle 6 between a direction (6*', 4>') within the cone, 
and the line of sight satisfies cosQ — cos9cos9' — sin9sin9' sin(f)' . Then the observed flux is 



where D = [F(l - pcosQ)]'^ is the Doppler factor, , /3 = (1 - F'^) V2^ = Dv' , is the 

specific intensity of synchrotron radiation at v' , and d is the distance of the burst source. 
Here the quantities with prime are measured in the comoving frame. For simplicity we have 
ignored the relative time delay of radiation from different parts of the cone. 

For the expanding jet, we have r = DVPcT oc DVT ~ 1), r' = D^cT oc DT, 
the magnetic field strength 5' oc F, the peak frequency of synchrotron radiation 



3. The emission from jet 
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— Du'^ oc DT^, and /'(j^^) oc n'^B'r' oc DV'^T. Assuming that the emission spectrum 
oc then r{v') = I'iy'Ji^y = I'i^'Ji^y^ j^i+a^^+sarp^^-a Therefore we 
have the flux 

F{u, 9) oc i/-"r2("-i)r3t,(^, r, a) (7) 



where 



g{e, r, a) = r d(j)' sine'de\i - pcosey 

Jo Jo 



(e+a) 



(8) 



In general, the value of g can only be calculated numerically. However here we consider the 
case 9j <C 1 and ^ -C 1, then cosQ ~ cos9cos9'. In this case we can calculate the value of 
g analytically under certain conditions. After complicated calculation we find g oc 
for r-^ < 9 < 9j, g (X ^|^-2(5+a)-2 ^j. p-i < ^ and ^ > ^ oc r2(5+«) for 9 < F'^ < 9j, 
and g oc r^'^^+^^+^^l for ^ < and < F"^. Therefore we have the results 



F{u, 9) oc < 



I/-°^2^-2(5+")-2r2(«-l)T3, for < ^, ^ > 



for ^ < < 9, 



(9) 



1-1 



^-a^2pl0+4aj.3^ for ^ < < T" 

For T <^ Th, the evolution is about F = Fo(r/rd)~^/^. Here we define the Lorentz 

factor, Tj = F(rj) = 9j^, then = (Fo^j)^/^rd, and the corresponding timescale 

Tj ~ rj/4F^c ~ 1.3(^)^-E'52^n]^^^^(^)^/^ day. Then the observed fiux at fixed frequency is 



F oc <^ 



rp-3a/2^ for T < Tj 



T-i"-i, for T,<T<T. 



(10) 



rp-2a-i^ for < T 

From above it seems that there should be two temporal index breaks in light curves. 
However, if we compare the values of Tj and Ti,, we will find that Tj <^ Ti,, i.e. the time 
interval between Tj and TJ, is very large, the beaming break is much earlier than the break 
due to sideways expansion. We know that in order to see the steepening of the light curve. 
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Tj or must be small, so in fact we can only see one temporal break. In addition, in the 
Rhoads' treatment, the effect of the sideways expansion on the F evolution was ignored 
when T < T^, however in fact, there is still some sideways expansion during this phase, 
so the evolution of F must be affected by sideways expansion when F ~ 9^^. Therefore, 
we expect that the evolution of F is continuous, and the transition from F oc T~^/^ to 
F oc T~^/^ is much smoother than previously claimed. 

In order to test our conclusion, we make a simple numerical calculation. We assume 
that the blast wave evolution is adiabatic, ignore cooling of the swcpt-up particles. We 
take the following initial parameters: Fq = 300, the electron distribution index p — 2.5. 
In Fig.l we present the light curve for different initial opening angle, the sohd, dotted 
and dashed hues represent the cases where = 0.1, 0.0174, 0.01 respectively. We show 
that when 6*0 = 0.1, the light curve is nearly not affected by sideways expansion, while 
when ^0 — 0.0174 and 0.01, the hght curves are steepened clearly, confirming our analytic 
conclusion. 

As for comparison, we also calculate the afterglow light curves of blast wave with no 
spreading. We know that, if the blast wave is highly radiative, the internal energy of the 
blast wave will be low, since it is converted to photons and radiated away, so the lateral 
expanding velocity will be very small, Cs « c, in this case, the sideways expansion is 
unimportantly small, and the blast wave can be regarded as no spreading. For highly 
radiative evolution, it has been shown that the Lorentz factor F oc T~^/^ (Wei & Lu 1998a). 
Taking the parameters as above, we calculate the afterglow light curves under this situation. 
Fig. 2 gives our results, the solid, dotted and dashed hues also correspond to the cases 
^0 = 0.1, 0.0174, 0.01 respectively. It is obvious that, when T <% (T^ ~ To(Fo^o)^/^), the 
light curves decay as a simple power law, and when T > Tc, the light curves deviate from 
the simple power law and show a steepening, and a clear break occurs at about Tg. Also we 
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note that only when Oq < 0.1, the steepening is obvious. Therefore we suggest that if we 
observe the sharp break in the GRB afterglow light curves, then it may indicate that the 
blast wave is highly radiative rather than adiabatic. 

4. Discussion and conclusions 

The GRB afterglows provide very good opportunity to study whether and how much 
the GRB ejecta are beamed. Rhoads (1997, 1999) has pointed out that the beamed outflows 
should diverge from the cone geometry and the sideways outflow of the shocked relativistic 
plasma would increase the front of the blast wave which leading to a fast deceleration. He 
also predicted that the afterglow hght curves should have a sharp break around T^. 

However, Moderski et al. (1999) have performed numerical calculation and shown that 
the break of the light curve is weaker and smoother than the prediction. Here we reanalyse 
the dynamical evolution of the jet blast wave, calculate the emission from the jet. Our 
calculations show that the main reason why the results of Moderski et al. being different 
from that of Rhoads is that the value of is very large when taking the parameters adopted 
by Modersh et al. (see eq. 5). Our formular (eq. 3) indicates that the evolution of Lorentz 
factor r with time T is continuous, changing the slope from -3/8 to -1/2 smoothly. In 
particular, if the value of is large, then the transition is much smoothly, in this case one 
expects that the sharp break will not exist. In order to test our analytic conclusion, we also 
make a simple numerical calculation, from Fig.l it is shown that only when < 0.1, we can 
observe the steepening of the light curve, which is consistent with our analytic conclusion. 

Our results are valid only if the remnant is still relativistic at time Ti,. Since at this 
time the Lorentz factor Ff, ~ |^^q \ this condition reduces to Oq < 0.1, i.e. the jet is very 
narrow. If > 0.1, then before the sideways expansion is important, the remnant has 
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already become non-relativistic. In fact, Dai and Lu (1999) have shown that even in the 
case of isotropic emission, the break in the hght curve can appear during the transition 
from ultrarelativistic to non-relativistic phase in the environment of dense material. 

It should be emphasized that in our calculation we have ignored the relative time delay 
of radiation from different parts of the cone, if this effect is considered, then the slope of 
the hght curve should be flatter than T'-(2a+i) (Moderski et al. 1999). However, all these 
calculation are based on the assumption that the material is uniformly distributed across 
the blast wave. In fact, it is more likely that the lateral outflow can lead to 9 dependent 
structure, with the density of swept material and the bulk Lorentz factor decreasing with 9, 
in this case the break in the light curve may become more prominent. 

We thank the referee for several important comments that improved this paper. This 
work is supported by the National Natural Science Foundation (19703003 and 19773007) 
and the National Chmbing Project on Fundamental Researches of China. 
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Figure Caption 



Fig.l The afterglow light curves for different initial opening angle, the solid, dotted 
and dashed lines represent the cases where = 0.1, 0.0174 and 0.01 respectively. 

Fig. 2 The afterglow light curves for highly radiative blast wave with no spreading. The 
solid, dotted and dashed lines correspond to the cases where — 0.1, 0.0174 and 0.01 
respectively. 
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